arXiv:1506.00657v2 [hep-th] 24Jul2015 


CPHT-RR017.0515 


Supergravity divergences, supersymmetry 
and automorphic forms 


Guillaume Bossard 1 and Axel Kleinschmidt 2,3 


1 Centre de Physique Theorique, Ecole Poly technique, CNRS 
91128 Palaiseau cedex, France 

2 Max-Planck-Institut fur Gravitationsphysik (Albert-Einstein-Institut) 
Am Muhlenberg 1, DE-lff76 Potsdam, Germany 

3 International Solvay Institutes 
ULB-Campus Plaine CP231, BE-1050 Brussels, Belgium 


We investigate the occurrence of divergences in maximal supergravity in various 
dimensions from the point of view of supersymmetry constraints on the U- 
duality invariant threshold functions defining the higher derivative couplings 
in the effective action. Our method makes use of tensorial differential equations 
that constrain these couplings. We study in detail the Fourier modes and wave- 
front sets of the associated automorphic functions and find that they are always 
associated to special nilpotent orbits. 
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1 Introduction 


The four-graviton scattering amplitude in type II string theory in D dimensions has an analytic 
part that possesses a low-energy expansion of the form 
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where £jj is the D-dimensional Planck length, s,t,u are the standard Mandelstam variables, and 
Cfc = (f|)/4) (s k +t k +u k ) dimensionless combinations of them. 1Z 4 denotes a specific contraction 
of the four polarisations and momenta of the four graviton states using the so-called tgfg tensor [1, 
2]. The most important objects in the above equation are the functions £(”,}) (Vo) that depend 
on the moduli Vd £ f?n_D/A’(Pn_D), where £n_c is the Cremmer—Julia symmetry group in 
D dimensions [3] and A'(f?n_o) its maximal compact subgroup. The functions £(^(Vd) must 
be invariant under the discrete U-duality group En_£>(Z) [4, 5]. Moreover, the functions have 
to satisfy differential constraints from supersymmetry [6, 7, 8, 9, 10, 11, 12, 13, 14] that can 
also be understood representation-theoretically [15, 16, 17]. The differential equations include 
Poisson-type equations of the form [9, 18] 
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for the first three functions £^o> and £ ( ( °o that are associated with R 4 , V 4 R 4 and V 6 i? 4 

type corrections to the supergravity action, respectively, and are of and A-BPS-type, 

respectively. 

Our interest here lies in the dimension-dependent source terms in (1.2c) for the V 6 i? 4 co¬ 
efficient £(q !) that we will determine by a new method in this paper. The source terms are 
related to perturbative divergences in supergravity. This is most clearly visible in the string 
perturbation expansion of the function £[qI) itself. The result given in [18] for the non-analytic 
part is 1 

£(o,l) = log 2 gs + y (| + 4?or) lo g58^ S D , 8 

20 5 

+ 5C(3) log g 6 S Dfi + —£ ( ( 0 5) 0) log g 5 SD ,5 + log g^D ,4 + ■■■■ (1-3) 

y 7r 

The term in So,a is related to a three-loop ultraviolet logarithmic divergence of the supergravity 
four-graviton scattering [19] and its value agrees with the field theory result of [20] as shown 
in [21, 18]. The terms in 60,5 and So ,4 are related to form factor divergences in supergravity. 
The precise coefficients in (1.2c) have only appeared recently in [18] and we will present here an 
independent derivation of these coefficients based on the tensorial differential equations of [10, 
11 , 12 ], 

There is a deep connection between constraints from supersymmetry and the Fourier modes of 
the functions <£’['” 1 [16, 17]. The Fourier modes arise in perturbative expansions of the type (1.3) 
where one exploits the periodicity of under discrete (Peccei-Quinn type) shift symme¬ 
tries of some axionic moduli. The non-zero Fourier modes in such expansions contain the 
non-perturbative corrections to the scattering process as they arise from instantons in string 
theory [5]. For the BPS-type correction terms in (1.2), only specific supersymmetric instantons 
can contribute and this puts restrictions on the structure of the Fourier expansion. Mathemat¬ 
ically, this is reflected in the so-called wave-front set of the Eu-oi^) invariant functions £{',?},)■ 
The wave-front set is a union of nilpotent orbits of the group Eu_o acting on its Lie algebra 
(see for example [22, 23, 24]). As Fourier modes can be associated with nilpotent orbits, the 
structure of the wave-front set captures the structure of the Fourier expansion and the correc¬ 
tion terms R 4 , V 4 R 4 and V 6 i? 4 can be associated with points on the Hasse diagram of nilpotent 
orbits [17, 10]. An important point that we will bring out in our discussion is that only so-called 
special nilpotent orbits are of relevance [17, 24], This will be discussed in detail for the case 
of 5*0(5,5) which is the Cremmer Julia group in D = 6 dimensions. We analyse carefully the 
Fourier expansion of certain Eisenstein series on this group that arise in the derivation of the 
logarithmic divergences in (1.2) and (1.3), presenting among other things the Fourier modes of 
the spinor Eisenstein series. 

This article is structured as follows. In section 2 we review Eisenstein series on symmetric 
space G/K as these are our main tools for constructing the correction terms We also 

introduce tensorial differential operators that are needed for writing the supersymmetry con¬ 
straints on the correction functions and Fourier modes in the subsequent chapters. In section 3, 

1 We have denoted the D-dimensional string coupling by go- When it is clear from the context which dimension 

we are in, we will often omit the subscript. 
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we present and solve the supersymmetry constraints in the case of Sffi) that is related to the 
three-loop divergence in D = 6 supergravity. Section 4 contains a new method for finding 
the divergent terms in other dimensions and derives the coefficients in (1.2c) and (1.3) from 
a particular ‘adjoint’ Eisenstein series on E^y Section 5 then analyses in detail the Fourier 
expansions of various Eisenstein series on 50(5,5) in connection to the supersymmetric cor¬ 
rections. Two supplementary appendices contain details on the adjoint E 8 (s) series and on the 
Fourier expansion of the spinor series of 50(5,5). 

2 Eisenstein series and tensorial differential operators 

In this paper, functions on the moduli space of string theory in D space-time dimensions play a 
central role. The moduli space is a symmetric space G/I\ with G = £n_o the Cremmer-Julia 
symmetry group in D-dimensional ungauged maximal supergravity and K its maximal compact 
subgroup. On this space, we will define Eisenstein series invariant under U-duality E\\_ d{Z) 
and tensorial differential operators that help to express the constraints from supersymmetry on 
functions on this space. 

2.1 Brief reminder of Eisenstein series 

We will use the following convention for Eisenstein series on a split real Lie group G, i.e. 
functions defined on the Riemannian symmetric space G/K, which are invariant under the 
arithmetic subgroup G(Z). Following the normalisation of Langlands [25], one can define an 
Eisenstein series for almost all (complex) weights A of G by the formula (see also [15, 26, 18] 
where also the other statements can be found) 

E( A,V)= e (A+p|H(7V)> (2.1) 

76 B(Z)\G(Z) 

where V € G is parametrised by the moduli of the theory in G/K. 2 The discrete subgroup G(Z) 
acts on G by left multiplication V i —> yV and the function H(V) picks out the logarithm of the 
Cartan torus part a in the Iwasawa decomposition V = nak of an element V E G. Therefore 
H(V) belongs to the Cartan subalgebra f) of G and can be paired canonically with the weight 
A € 1}*. The element p E fi* in (2.1) denotes the Weyl vector (half the sum of the positive roots). 
The series is absolutely convergent for 

(Re(A)|a) > (p\a) for all a > 0 , (2-2) 

for all positive roots, and extends to a meromorphic function of the weight A over !)* [25]. The 
function E( A, V) is G(Z) invariant and satisfies the Laplace eigenvalue equation 

A£(A,V) = |«A|A> - (p\p))E(X,V). (2.3) 

2 The Eisenstein series is by construction spherical meaning that E(X,Vk) = E( A,V) for all k £ K and can 
therefore be viewed as a function of the moduli in G/K. 
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in terms of the standard bilinear form (•]■) on I)* which is normalised such that long roots a 
have length {a\a) = 2 . 

For the particular series arising in this work it will be convenient to parametrise A as 

A = 2u j-p. (2.4) 

The advantage of this notation is that for maximal parabolic Eisenstein series the weight w G t)* 
thus defined is proportional to a fundamental weight of G. More precisely, we expand u on the 
basis of fundamental weights Aj (i = 1,..., rank(GQ) as 

dim lj* 

UJ = ^2 S i A i- ( 2 - 5 ) 

i =1 

A maximal parabolic Eisenstein series then has only one non-zero Sj. Using (2.5) we denote the 
Eisenstein E( A, V) alternatively by a labelled Dynkin diagram using the numbering conventions 
of Bourbaki (identical to those of the LiE program [27]). For example, for 50(5,5) of Cartan 
type -D 5 we will write 


E 


( 2 . 6 ) 


and will always suppress the dependence on the coset representative V £ G. An example of a 
maximal parabolic Eisenstein series for E 8 ^ then would be 


E r o l- (2-7) 

[s 0 0 0 0 0 oj 

This is defined for almost all complex s (by analytic continuation). The value s = | corresponds 
to the R 4 function (1.2a) and s = | to the D 4 R 4 correction (1.2b) [28, 15]. 

We will also sometimes refer to a maximal parabolic Eisenstein series by the representation 
the relevant fundamental weight A $ corresponds to. In this terminology 


E 


will be called a vector Eisenstein series of 50(5, 5) (for any s) and 


( 2 . 8 ) 


E r o 1 (2.9) 

[o 0 0 0 0 0 sj 

an adjoint Eisenstein series of E s r 8 ) (for any s). 

We must warn the reader that in this paper we will always consider the Eisenstein series in 
the Langlands normalisation (2.1), whereas one often finds the lattice sum normalisation that 
differs by a factor of 2 £( 2 s) in the literature. 
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2.2 Functional relations and constant term formulas 

Eisenstein series satisfy almost everywhere the functional relation 


E( A,V) = M(w,X)E(wX,V), 


( 2 . 10 ) 


where w is an element of the Weyl group W = W(G') and the intertwining coefficient (sometimes 
also called reflection coefficient) is given by 


M(w, X) 


n 

a>0 

wa<0 


e((A| a)) 
a(\\ a) + iy 


( 2 . 11 ) 


where the product is over all positive roots a that are mapped to negative roots by the Weyl 
word w. The completed Riemann zeta function 


£(s) = 7T -S / 2 r(s/2)e(s) 


( 2 . 12 ) 


has simple poles at s = 0 and s = 1 with residues —1 and +1, respectively, and vanishes nowhere 
on the real line. It satisfies the functional relation £(s) = £(1 — s). 

Of use to us will also be constant term, formulas that express the integration over some of the 
variables in the unipotent part U of a parabolic subgroup P = LU C G in terms of automorphic 
functions on the Levi part L. In physical terms, the constant term formula expresses the result 
of averaging out certain axionic moduli and thus projecting to the zero-instanton charge sector 
for instantons charged under these axions. A parabolic subgroup is the product of an abelian 
subgroup GL( l) xn and the semi-simple component of the Levi subgroup. The constant term 
formula projects onto the perturbative part in the moduli parametrizing this abelian subgroup 
GL( l) xn , which are generically combinations of the string coupling constant and the radii of the 
compactification torus. We will label maximal parabolic subgroups (corresponding to n = 1) 
by P, = LJJi, where i denotes the node of the Dynkin diagram of G that has to be removed 
to obtain the Dynkin diagram of the Levi subgroup Li (more precisely that of the semi-simple 
part Gi := [Li, Li] of Li = GL{ 1) x [Li, Li] where [Li, Li] denotes the commutator subgroup). 
For a maximal parabolic subgroup Pi the constant formula can be written as [29, 15, 26] 

fE(X,uV)du= M(w,\)e {{wX+p) Wi mv)) E Gi ((wX) ±i ,Vi). (2.13) 

Ju i WGWiXW 

Let us explain the notation in this formula. The integration domain is Uj := Ui(h)\Ui = 
(G(Z) n Uj.)\Ui which is the fundamental domain of the discrete shifts in the unipotent group 
C/j(Il) to restrict the integration to a single period. In the simplest case of a one-dimensional 
unipotent and in a convenient normalisation one has U} = Z\R = [0, 1), the unit interval giving 
rise to the notational superscript 1 . Since we are averaging over the Ui dependence, the result of 
the integral can only depend on the variables parametrising Li = GL( 1) x Gi. The dependence 
on the two factors is separated on the right-hand side, where the dependence on V* € Gi is via 
an Eisenstein series on the group Gi and the dependence on the GL( 1) Cartan torus factor is 
written in terms of the exponential prefactor. A given weight A of G can be decomposed into 
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a component parallel to the fundamental weight A* (by orthogonal projection) and remaining 
components orthogonal to it: 

A«Ay. + A ±i , Ay ( = ^ Ai | A .^ A * • ( 2 - 14 ) 

When Ay. is contracted with an element of the Cartan torus it picks out only the component 
along the GL(1) factor in the Levi subgroup L* and therefore the exponential prefactor stands 
for some power of a variable on GL( 1) whose normalisation we will choose to give it an easy 
physical interpretation. The component Aj_ ; is then a combination only of the simple roots of 
the subgroup Gi C G and can therefore be used to define an Eisenstein series on the group Gj. 
The sum in (2.13) is over the quotient of the Weyl groups of G and G, and the (numerical) 
coefficient M(w, A) is precisely the intertwiner defined in (2.11) and hence given by a quotient 
of completed Riemann zeta functions (2.12). In keeping with our notation we will typically 
suppress the moduli dependence in the Eisenstein series and label the weight A in terms of its 
Dynkin diagram representation as in (2.6). (There exists a different constant term formula when 
U is the maximal unipotent N [25] but we will not need it here.) 

As an example for (2.13), we can consider the following constant term integral for the vector 
Eisenstein series of 50(5,5) 



= 2s 2B - 8 g(2s-4)g(2s-7) 

9 ^ 9 £(2 S )£(2 S -3) 


+ g 


-i 


g(2* ~ 1) F 

£(2 s) 


(2.15) 


where g denotes a coordinate on the GL(l)-part of the Levi subgroup L\ = GL(1) x 50(4,4) 
of the maximal parabolic P\ associated with the first (left-most) node of the 50(5,5) diagram. 
50(5,5) is the Cremmer—Julia group in D = 6 space-time dimensions and the constant term 
formula above corresponds to a string perturbative expansion since it preserves the (chirality 
preserving) T-duality group 50(4,4, Z) in D = 6 and that is why we have labelled the expansion 
parameter by the six-dimensional string coupling g = g g. The labelling of the 50(4, 4) Dynkin 
diagram is chosen such that the first node corresponds to the vector representation of the T- 
duality group, noting that the RR moduli define a Weyl spinor of negative chirality of Spin( 4,4). 
As is well-known [30, 28, 15], the above Eisenstein series is related to the six-dimensional R 4 
correction by 


5 ( ( 0 6) 0) = 2C(3)^ foo o| • (2.16) 

For the particular value s = §, the middle term in the constant term formula (2.15) vanishes 
due to the properties of the completed Riemann zeta function and one obtains therefore 

5 (0,0) = 2C(%' 3 + 4C(2 )g~ 1 E^ q ; j , (2.17) 

corresponding to the correct tree-level and one-loop contributions to the |-BPS coupling R 4 . 3 
Manipulations of this kind will be central for evaluating the supergravity divergences in various 
dimensions that arise from poles in Eisenstein series. 

3 As usual, the symmetry is made manifest in Einstein frame, explaining why the powers on the string coupling 

are shifted from the string frame values g~ 2 and g° for tree level and one-loop. 
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2.3 Differential equations 


Eisenstein series are eigenfunctions of all Casimir differential operators on G/K for almost all 
values of the weight A € f}*. To define the differential operators on G/K it is convenient 
to take a specific representation of the coset representative V € G in terms of coordinates 
<h, as for example the one, = (<)>, er), associated to the Iwasawa decomposition, such that 
V(<h) = V(4>, ct) = n{a)a{cj))k. One can define in this way the symmetric space connection Q and 
its vielbeins P from the components of the Maurer-Cartan form restricted to the Lie algebra t 
of K and its orthogonal complement in g, i.e. 

V($) _ 1 dV($) = Q^)d^ + P^)d^ , Q e 6 , Pegef. (2.18) 

The group K defines the structure group of the symmetric space, and one can modify the 
reference frame by arbitrary functions fc(<h) such that 

Q -> + L($)” 1 Qfc(^) , P -»■ fc($) -1 .PIfe($) . (2.19) 

The Riemannian metric on the symmetric space G/K is defined for some appropriately nor¬ 
malised G-invariant bilinear form 


G^) = (P^),P,m , ( 2 - 20 ) 

and permits to define the inverse vielbeins through its inverse 

V 11 = G lw P,j . (2.21) 

One defines the covariant derivative in tangent frame V, as the differential operator acting on 
any tensor function //? K ( < h) on G/K in an arbitrary representation Rk of K and transforming as 
fn K —> 'KR K {k~ 1 )fpt, K . The differential operator takes values in the R-represent at ion (g©t)(g )Rk 
and is defined by 

P/ Rx (<h) = © (9 M + 7T R K {Qv))fR K (*) , (2.22) 

where ttr k (X) is the Lie algebra element X e f in the representation Rk acting on fn K (T). 
In particular, for a function defined on G/K such as £(//,/, one defines recursively the n th order 
differential operators 

T>m) = V^W) 

f v d u tt$j) + [Q tt ,V v ]d v tt$j) 

V®V®Vf{§) = V*®dp(v®Vf($)}+V»®[Qr,V v ]®(d l ,Vf($) + [Q v ,Vf($)]) 

+V» © V v © [Q„, {d v Vf(< D) + [Q„, P/($>)])] (2.23) 

which is valued in the n th tensor power of the Lie algebra component g © t. 

The differential operators can be written in an arbitrary representation R of g by writing the 
coset element in the representation R. Then powers of the differential operator are mapped 
to powers in the representation R, such that one projects these differential operators valued 
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in the tensor algebra to the enveloping algebra of g associated to this representation, and one 
writes then 


Dfl/($) = 7T R (Vndjm 

D i/(<h) = +[n R (Q^,V R fm) 

D Z +1 W) = 7T R (V^[d,{n^m) + [7r R (Q,),T>Sfm) , (2.24) 

which defines a matrix of differential operators in an explicit matrix representation R. Doing so 
one considers by construction the restriction of these differential operators to specific irreducible 
representations of I \, with the technical advantage that it becomes relatively simple to compute 
the explicit form of these differential operators in specific parabolic decompositions. 

A generic character satisfies by construction tensorial differential equations in 

some irreducible representations of K , depending polynomially on the weight vector u. For uo 
such that the bound (2.2) is satisfied, the associated Eisenstein series is absolutely convergent, 
and it follows that it also satisfies the same tensorial differential equations for almost all u 
by analytic continuation. For some sub-classes of weight vectors ui, in particular when the 
latter is proportional to a fundamental weight, i.e. uo = sA t . the generating character satisfies 
stronger differential equations that can often be rewritten as characteristic equations in a given 
representation R, i.e. 

V u , R (B R )e 2 ^ H ^ = 0 , (2.25) 

for a polynomial r(Dr) in the covariant derivative D# that depends polynomially on the 
weight vector oj (see [12] for some examples). Whenever this sub-class defines a domain intersect¬ 
ing with the domain of convergence of the corresponding Eisenstein series E{ A, V) at A = 2a; — p 
(2.2), it follows that the latter satisfies 

V u ,r<Pr)E(\,V) = 0 . (2.26) 

Because E( A, V) is analytic in u for almost all c v and the differential operator V U) r{T>r) is 
analytic in <n, this equation is then satisfied in general for uj. 

It may also happen that a generating character satisfies a stronger characteris¬ 

tic equation for an isolated u, for which the corresponding Eisenstein series is not absolutely 
convergent. Then one cannot directly conclude that the Eisenstein series satisfies itself this 
stronger characteristic equation. In particular, when the differential constraint on the character 
is associated to a nilpotent orbit that is special in the sense of [31] the corresponding Eisenstein 
series will also satisfy a weaker constraint associated with a larger special orbit [32, 24]. 

Let us for this purpose consider the example of E 7 ^ 7 ) Eisenstein series. It was computed in 
[12] that 

D^e 2 ' 5 < A 'l"(V)) = ^(2 a -17) +6 ^ D56e2a(Al | H ( V )> > 

DiV s(A7|H(V)> = a(a - 9)D 133 e 2 *< A ^W> , (2.27) 
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and because the corresponding Eisenstein series are respectively absolutely convergent for s > hy 
and s > 9, one concludes that for almost all s ( i.e. away from the poles) 


D 5 6 ^ 


D 133 -E 


s 0 0 0 0 0 

0 

0 0 0 0 0 s 

4 d4 


s(2s - 17) 


+ 6 ) D 56 -E 


56-^1 0 

sO ocool 


= s(s — 9)Di 33 £r 0 n 

0 0 0 0 0 s 


(2.28a) 

(2.28b) 


In type II string theory, the V 4 £ 4 threshold function is conjectured to be £(5)£[s 0 ° 0 0 0 | [15], 
and must satisfy by supersymmetry [ 10 , 11 ] the two equations 


d 5 3 6 £ 


0 

#00000 


D133-E7 _ 0 


#00000 


= —20D133E 


(2.29) 


The first equation is obviously satisfied by C(5)£[s 0 ° 0 0 0 J. For the second equation, one has to 
use that, because of the Langlands functional identity (2.10) 


C(5)£r 0 

0 0 0 0 0 


8C(8) 

157T 


E \ o I , 

[0 0 0 0 0 4j 


(2.30) 


this function is a special case of both the adjoint and the fundamental Eisenstein series. From 
equation (2.28b) for the fundamental Eisenstein series one then sees that both supersymmetry 
constraints are fulfilled. Moreover, one can in this way understand that this function admits a 
wave-front set associated to the next-to-minimal nilpotent orbit of £7 [17]. 

By supersymmetry, the type II string theory R 4 threshold function must satisfy the stronger 
differential equation [ 10 ] 

Dt 1 C) = ( 2 - 31 ) 


and the conjectured solution £$ 0 ) = 2C(3)£[a 00000 ] [28, 15] indeed solves this constraint [10]. 

However, one must in general be careful when the Eisenstein series is outside the domain 
of absolute convergence. We will see in section 5.2 that even though the character e 3 (A 2 |-ff(V)) 
of the adjoint series of 50(5, 5) satisfies a certain stronger constraint the associated Eisen¬ 
stein series (and its Fourier coefficients) do not. This also happens for the adjoint series of 
£ 7 ( 7 ) and the character e 8 (Ai|.ff(V)) that satisfies an additional quartic differential equation in 
the [2,0, 0,0, 0,0, 2] of 51/(8), while the corresponding adjoint function at s = 4 violates this 
constraint. In appendix A we show the same for the adjoint £s( 8 ) series. 


3 Supersymmetry constraints on £ (01) 

It was shown in [12] that the V 6 £ 4 threshold function £^o,l) decomposes into the sum of two 
distinct functions associated to two different supersymmetry invariants for D > 3. The two 
invariants are distinguished by higher point R-symmetry violating couplings. One function sat¬ 
isfies a homogeneous differential equation and is conjectured to be an Eisenstein series, whereas 
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the other satisfies an inhomogeneous equation as was first argued in [ 8 ]. For instance, in six di¬ 
mensions, Erj = SO( 5, 5) and K = 50(5) x 50(5), and the threshold function £( 0 ,i) decomposes 
as 

c . 4C(6)«8) p 

£»,.) - A0,1) + 2^(4) *[„ „ „ • ] • 

(The hat here denotes a regularised spinor Eisenstein series that will be defined in (3.6) below.) 
The supersymmetry analysis only constrains the function appearing in the Wilsonian action, 
and one must consider possible anomalous corrections to the corresponding differential equation 
whenever there are logarithmic divergences in the theory. It turns out in particular that the 
four-graviton amplitude diverges at 3-loop in six dimensions [20], and supersymmetry therefore 
only constrains the function Jy 0 ,i) to satisfy 

2 70 


A-F'(o,i) = - (2C(3)^ 3 00 0 ]) + y crC(3) , 


where c\ is a constant yet to be determined, and [ 12 ] 

( 0 , 1 ) = —Dle-^";"( 0 , 1 ) — 2£(3) E 


3 0 0 d 


D? n ^ 


10 -Ho.i) — ^^lO-Tyo,!) ~ 2C(3) E | 2 0 0 0 


1D16-E 


DinE 




'““boo 8 


(3.2) 

(3.3a) 

(3.3b) 


Here, Di 6 refers to the covariant derivative valued in the (chiral) spinor representation, and 
Dio t° the covariant derivative valued in the vector representation, according to the notation 
introduced in the preceding section. One can for example write Dig and Dio as explicit matrices 
of differential operators as follows 


D 


V 


10 — 


ab 


1 


a.b 


0 


D 16 = 2 V abT7 


(3.4) 


where T> ab is the covariant derivative as a (5,5) tensor of 50(5) x 50(5), and we label the 
vector indices of the first factor by a (ranging from 1 to 5) and those of the second factor by a. 
7 °, 7 “ are the Spin( 5, 5) gamma matrices in a (fixed) Majorana-Weyl representation and 50(5) 
indices are raised and lowered with the flat metric. 

It is important to note that equations (3.3) are invariant with respect to the exchange of 
chirality, and read in 50(5) x 50(5) covariant notations 


"abode 


V a a V b b V c c E 0tl) = 0 


v a£ v dc v di E {0A) = -V ab E (0A) - 2C(3 ye 


io 0 0 0 ] V ab E \ ioo l 


On the other hand, one computes using 1 


E 


= lim (E 


m 1 


e-^oV [°°° 4 + e J £(6)£(8) 2e 


(3.5a) 

(3.5b) 

(3.6) 


4 The regularised function is only defined up to an arbitrary additional constant, which is associated to the 
ambiguity in defining the separation between the local and the non-local components of the effective action in the 
presence of logarithm terms. 
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that [12] 


AEr 


[° 0 0 ° ] 5 ae)e(8) ’ 


e(3) 


E Uoo o r -± £abcde£ ^ vdp e* o o 


which is not invariant under the parity transformation in the T-duality group. 


(3.7a) 

(3.7b) 


3.1 Solution using string perturbation theory 

We will now construct the solution to the tensorial differential equation (3.3) using string per¬ 
turbation theory. At the end of the derivation we will argue that the differential equation (3.3) 
does not admit cusp form solutions, and the method provides the full non-perturbative solution. 
According to string perturbation theory, the function Jqo.i) decomposes as 

•F(o,i) = e- 6< ^ £ c 2 </> + 0(e- e ' 0 ) , (3.8) 


where = g is the string theory effective coupling constant in six dimensions, and the additional 
term linear in <f> must be added to take into account the non-analyticity of the threshold function 
due to the 3-loop divergence [19]. To solve these equations in the string perturbation theory 
limit, we need the explicit decomposition of these differential operators acting on a function T 
of the dilaton <f> and the SO( 4,4) scalar fields' 5 


Dio-F = 


0 

0 


0 

Ds a 

0 


0 

0 

~lp<t 


T 


~Di§E = 


4W 


ich, + Dg 


0 


— \d(j> + Dg 


T 


(3.9) 

Note that we use the embedding of 50(4,4) as the T-duality group in string theory, referring to 
the property that the 16 vector fields in six dimensions are associated to 8 NS and 8 RR fields 
(respectively in the 8 and the 8 C of Spin( 4,4)). One computes that the second order differential 
operator defined as in (2.24) is 


= 


Di= 



0 

D« 2 + id. 


I u <t> 


T 

J J 


0 \d£ + 

+ \®<t> + \ {®<t> + l) D 8 + D 8 2 


0 


+ \d<t> “1(^0 + 1 ) D 8c + Dg 


T , (3.10) 


5 Here \ d$ + D 8 is understood to be + Ds- We will never write explicitly the identity matrices in the 

following. 


11 





and finally 


Dio-F 


jy 3 tt 


( H 3 + 0 0 2 + f^-iA 0 0 \ 

0 D 8 3 + i^D 8a 0 Lf , (3.11) 

V 0 0 -H 3 - 9 0 2 -§^ + ±A ] 


A/) 3 I In 2 _ 

64°V + 8 °<t> 

D — —r) 3 
U 64 °<f> 


§ A + il6^^ + T^<t> + l) D s + f {d<f> + 2)Dg 2 + D g 3 0 

“ + l A + (il d 0 + TT^ + l) D 8c — f ip<!> + 2 ) D 8 2 c + D 8c 


T . 


We recall that the square notation is a short-hand notation for the definition (2.24), explaining 
the additional lower order differential operator contributions associated with the terms involving 
the connection. By construction 

a £> 5 -F = trDjoJ 7 = -trD, 2 6 F = ^An 4 + -d£ + 4 . (3-12) 


To solve these differential equations, we will make use of the particular solutions 


B i E [ s o g ] 

- { s ( s - 3) + |)D 8 Fj s q o j . 

s(s - 3) p 
~ 4 !-.;]• 

Dc? E\ o i 

_ s ( s “ 3) p 

4 [ 00 °] ’ 


D ^[-oS] 

Do 3 E\ o i 
Oc 00 a 

= ( s ( s - 3) + |)D 8c Fj o Q o j , 

(3.13) 

D «4o:] 

»(»- 3) p 
“ 4 

D 8 ^[ooO] 

_ s ( s ~ 3) p 

4 [ 00 °] ’ 


and equivalently for E ^ q s j. 





We will now analyse the differential equations order by order in string perturbation theory, 
i.e. for the various Using equation (2.17) with g = e^, one can now compute 


3 3^ \2C(3) 2 


Dio - ~D 


2 10 7 3 


6 <t> 


r, 3 3 n 2C(3) 2 w 
° 16 — 4 Dl6 / 3 6 


such that 


Similarly, one computes that 


/ 1 

0 

° ^ 


3C(3) 2 0 

0 

0 

e~W = -2C(3) 2 e“ 3 ^D 10 e" 3 ^ , 

V ° 

0 

- 1 J 


5« 3 ) 2 ( 0 

0 

-l 

J e~ 6 * = -2C(3) 2 e" 3</, Di 6 e" 30 , (3.14) 

tt(°) _ 

^(0,1) 

2C(3) 2 

3 


(3.15) 



V,) 4C(2K(3) e-^ 

2 / 3 





4C(2)C(3) ^ 

3 



4C(2)C(3) 


/ 2 0 

0 —Dg 

^ 0 0 


0 \ 

0 

-2 


E 


e in r o 
10 o 


-4C(2)C(3)D 10 e- 4 ^ io oj , 
4C(2)C(3) 

-4C(2)C(3)Di 6 e- 4 ^ 


1-D 8 0 

0 — 1 — Dg 


-W E 


o > 


(3.16) 
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such that 


irC 1 ) 

^(0,1) 


4C(2)C(3) 


E 


(3.17) 


Note moreover that there is no homogeneous solution to these differential equations with the 
corresponding factor of the dilaton, such that these solutions are unique at these orders. The 
2-loop contribution satisfies the more complicated equations 


0 


Di 3 o-TDifl)e-% = 


(o.i) 


3 _ I A 
2 4^ 

0 Do 3 - 2D 8 


0 \ 


0 0 -| + |A J 


-2(j> tt(2) 
e (0,1) 


and 


-D 10 (2C(2) e -^ [io „|) ; 


2 


1 0 0 \ 
0 -D 8a 0 

0 0-1 J 


(2C(2)e-^ io oj) 2 (3.18) 


0 


3_1 A + D 3-2D 8 

— i + |A + D 8 3 — 2D 8 


0 


'8 C ~ Z1J 8 C 


-20 x-(2) 
e -'(0,1) 


= -D 16 (2C(2)e-^j io 

from which one deduces that 

a tt( 2 ) _ c 7r( 2 ) _ 

(0,1) — U -' (0,1) 


1 


Ds 


-1-D 8 


(2C(2)e-^ io oj) 2 , (3.19) 


2 

, 0 I ) J 


(Ds- -2D 8i )^ = -D 8i (2C(2)^ io oj) 2 , (3.20) 


where 8^ stands for the three fundamental representations of Spin( 4,4), the vector and the 
two Weyl spinor representations. It is important to note that this tensorial equation is triality 
invariant since 



= E 



(3.21) 


This property is crucial in the conjecture proposed in [18] that the function J 7 ^, i) is triality related 
to the genus two V 6 i? 4 threshold function in five dimensions. The differential equation for 
admits an 50(4,4, Z) invariant homogeneous solution, such that the differential equation only 
determines the correct function up to 


7r( 2 ) 

^(o,D 


_T*( 2 ) 

(0,1) part 


l 87r 6 C(5) 2 
496125C(7) 3 [oa ° 


(3.22) 
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Let us now derive the 3-loop contribution. One cannot disentangle the cj) independent function 
from the logarithm term linear in 4> hr the differential equation, so we consider 


such that 


D,t, - -D 10 

5C£3) _ 1 A 

2 02 4 Z - X 

0 


TT (3 ) , 

■> (0,1) ' 


10C(3) 


-c 2 


0 


0 

0 


D 8 3 a - §D 8 


5C(3) 


0 

0 

c 2 + jA 


f (3) — 0 
(0,1) ~~ u 


Die - -D 


16 


A) + 


10C(3) 


-c 2 


5C(3) 


c 2 - |A + |Dg 2 + Dg 
0 


5C(3) 

8 


c 2 + — §D g 2 + D 


A) = 0 , (3.23) 


D 8 3 a ^ 0 3) i) = xD s . T, 


(3) -p> 2 -77(3) 5C(3), ■n 2 tt( 3) _ hC(3) 

D 8 -0(0,1) a c 2 > L' 8c *' (o,i) A C 2 


(o, 1 ) > —a~ io,i) ^ “•s > “o 0 - (u,i) ^ 

These differential equations are solved by the regularised Eisenstein series 

£(3) 1 


E r 3 1 = liml Et 3 , «i _ . _ . . . 

■ 0 0 J 0 0 + J f( 4 )f( 6 ) 2 e y 


(3.24) 


(3.25) 


which by construction (3.13) satisfies 


d 8 3 £ 




d* 2 e 


3 g(3) 
8£(4)£(6) 


D b 2 A 0 3 


3 g(3) 

8£(4)£(6) • 


(3.26) 


Assuming that there is no cusp form satisfying these differential equations, one can argue that 
this is the unique 50(4,4, Z) solution with a sufficiently fast fall-off at the boundary of moduli 
space. To show this we use the property that the differential equations determine the eigen¬ 
values of all the Casimir operators. It follows that the general solution on the maximal torus 
(i.e. the infinitesimal quasi-character fixed by the weight A in (2.1)) is uniquely determined up 
to Weyl reflections, and so is the general Eisenstein series solution. The spectral decomposition 
of automorphic forms is into a continuous part (given by Eisenstein series) and a discrete part 
(corresponding to cusp forms and residues of Eisenstein series) [33] . One can verify that there 
is no residual spectrum for the weight A defined by the above tensorial differential equations. 
Therefore the assumption on the absence of cusp forms implies that the 50(4,4, Z) Eisenstein 
series solution is the unique solution. The existence of a cusp form satisfying these differen¬ 
tial equations would not affect our conclusions as the cusp form would not contribute to the 
perturbative series on which we base our analysis. 

As we will discuss in the following, comparison with higher dimensional string theory com¬ 
putations allow to determine J 7 ®!) to be 


X*(3) _ 

-Mo,!) “ 



(3.27) 
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such that C 2 = 1. Using this function, one computes that 

Af 


3 1 

+ 10C(3) A 

V 27 [ 00 oJ 

3 J 


70C(3) 


(3.28) 


and therefore c\ = 1 . 

By arguments similar to above, there is no solution to the differential equation at higher 
order in string perturbation theory. This establishes the expected non-renormalisation theorem 
that the function £( 0 ,i) is exact at 3-loop in perturbation theory. Combining these results one 
deduces that 

.Ann\rii\ 16C(8) 


_ 2C(3 ) 2 6 , w 4C(2)C(3) 

3 e + e 3 


E 


+ 


-E 


+ 


4C(6) 

27 


E 


189 
H +E 


i p - 2</> tt( 2 ) 

' e •> ( 0 , 1 ) 


+ 5C(3)0 + O(e- 


(3.29) 


in perfect agreement with [18]. 

We would like now to argue that the differential equation (3.3) does not admit cusp form 
solutions, and therefore determines the complete non-perturbative function uniquely. Cusp forms 
are square integrable functions £ on the modular domain Eg — G(Z)\G/K and eigenfunctions of 
the Laplace operator. Due to their cuspidal nature they have strictly negative Laplace eigenvalue 
since 


d[x£A£ = — dfj, |V £| 2 < 0 . 


(3.30) 


JTg j Fg 

Therefore, there can be no cusp form solution to equation (3.3). 

Let us summarise what we have achieved through this computation. Consistency with 
higher dimensional results permits to determine that the 3-loop contribution Eff^ is equal 
to ^ 27 *) E [o o g J + C 2 10< ^ 3 ) 0 , and we have shown that the tensorial differential equation (3.3) is 
strong enough to determine C 2 and the anomalous source term in the Laplace equation, i.e. c\. 


3.2 Alternative derivation of non-analytic terms 

We will now argue that one can alternatively use the regularised Eisenstein series E\o 5 0 ° J to 
probe these properties, although the latter does not appear in the complete threshold function 
E{q \j . The adjoint Eisenstein series satisfies in general that 

2s(2s - 7) +3 


D AE 


A E 


-Di6 E 


Dio E 


s(2s- 7)+ 3 
01 —---UlO^ 


= 2s(2.s-7)E [os0 o 


0 J 

(3.31) 


050 ° is a ho- 


such that the regularised Eisenstein series with the normalisation 4 27 gf 2 ) 6 * gfi) 8 ' > ^ 
mogeneous solution to equation (3.3) that reproduces precisely the anomalous term in (3.2). 
Therefore, the ‘regular’ function 6 

4£(5)C(6)£(8) 


• 77 K — 77 _ 

•r ( 0 , 1 ) — -Do.i) 


27^(2)C(4) 


-E 


, 0 > 


(3.32) 


6 We call this function regular because the absence of anomalous term on the right-hand side of (3.33) suggests 

that it lies in a continuous family of functions that would be regular at this point. 
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defines a particular solution to (3.3) that is a solution to the Laplace equation 

AJ-* 1) = -(2C(3)^3 0o oj) 2 • (3.33) 

By construction, the 3-loop contribution to E(o,\) is the corresponding constant term of the 
regularised Eisenstein series E\o 5 0 ° J, as well as the non-analytic term linear in cj) associated 
to the logarithmic ultra-violet divergence. So instead of computing explicitly the solution to 
the tensorial differential equation (3.3), one could have instead used the property that the 
regularised Eisenstein series E\o 5 0 ° J is the unique automorphic homogeneous solution to this 
differential equation, and use the Langlands constant term formula to derive the solution (3.28). 
This computation can easily be generalised to arbitrary dimensions, and this is the approach we 
shall expand on and follow in the next section to derive the non-analytic components of £( 0 , 1 ) in 
dimensions four and five. 


4 Divergent pieces across various dimensions 


Our strategy for determining the coefficients in (1.2c) outlined at the end of the preceding section 
can also be stated as follows. We decompose the complete threshold function as 


c( D ) _ -WD) R . nj(D) 

C (0,l) _ J (0,1) rL (0,l)> 


(4.1) 


where 


is a particular (regular) automorphic solution of the non-anomalous inhomogeneous 


equation 



6(14 — D)(D — 6 ) 
D - 2 


-r-(U)R _ 

•Oo.p — 



(4.2) 


and the tensor equation associated to the adjoint Eisenstein series with s = 6 , | and ^ in 

dimension D = 3, 4, 5 and 6 , respectively [12]. It is very important that we choose 
such that it does not include any three-loop contribution in string perturbation theory. We can 
always do this because the three-loop contribution is generally the solution to a homogeneous 
equation. In fact the anomalous contribution is always associated to the three-loop contribution: 
in six dimensions through the 3-loop ultraviolet divergence, in five dimension through the 2- 
loop divergence of the 1 -loop R 4 type form factor, and in four dimensions through the 1 -loop 
divergence of the 2-loop V 4 !? 4 type form factor. 

From this definition it follows that is a solution to the anomalous (almost) homogeneous 
equation 



6(14 — D){D — 6 ) 


D -2 


n 


(O) _ 
( 0 . 1 ) — 


^ D ^Ad,s= (6-uK 1+ g) 1 


(4.3) 


where E\d,s is the adjoint Eisenstein series of E7 ii_d for D < 6 (which is 1 at D = 6 ), and 
/3d is the corresponding numerical coefficient, which will be determined in the sequel. Such a 
right-hand side must be considered whenever power-counting a priori allows for a form factor 
or a genuine amplitude logarithmic divergence in supergravity and the occurrence of the adjoint 
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Eisenstein series is fixed by the differential equations of [12]. Supersymmetry Ward identities 
moreover imply that such a correction can only occur if the two functions satisfy to compatible 
differential equations and in particular if the Laplace eigenvalues are identical. Investigation 
gives that this only occurs for D = 6, 5, 4. is the sum of two respective solutions to the 

two tensorial differential equations associated to the two distinct V 6 i? 4 * type invariants [12]. Up 
to cusp forms, the automorphic solution (with appropriate fall off at the boundary of moduli 
space) to these differential equations is unique, and the T-duality invariant solution is also 
unique at a given order in string perturbation theory. This implies that the general T-duality 
invariant solution to the differential equations relevant at 3-loop is necessarily the string theory 
limit of the fully automorphic solution. We can therefore assume without loss of generality that 
is the sum of two regularised Eisenstein series solving the two relevant tensorial differential 
equations, i.e. an adjoint Einsenstein series with s = , 6, | and | in dimension 3, 4, 5 and 6, 

respectively and a series in the fundamental representation associated to the last Dynkin node 
in the Eu_d convention with s = 14 ^ D , which define the same function in three dimensions. 

Note that in general the adjoint Eisenstein series admits an expansion in the perturbative 
string theory limit that is incompatible with string perturbation theory. There is for example 
always a term that would formally contribute to a — ^-loop correction. One can easily understand 
this in type IIB supergravity in ten dimensions, where one can define 11 such that it does 
not include a 3-loop contribution, in which case 


n 


( 10 ) _ 
(o,i) — 


4C(6)g(8) 

27£(7) 


E[ 4] (fi) 


4C(6)g(8) 

27£(7) 


^2 + 




(4.4) 


where the first factor is a spurious — ^-loop contribution. These spurious contributions all occur 
at negative order, and one understands that they are precisely compensated by the particular 
solution 


4.1 Three-loop divergence in D = 6 

The following (formal) Eisenstein series are the two automorphic homogeneous solutions to the 
two respective tensorial differential equations (3.5b) and (3.7b) 


E 


and E 




(4.5) 


Both of them are singular and need to be regularised.' A combination that provides a regular 
limit is 


T-lfoi) = ° li m I Er o 1 + TTxf^r 

(0,) [°°o 4 - £ J f(2) [o 


£(5) 


y+ e o 


(4.6) 


with a some constant that will be determined in the sequel. The relative coefficient here was 
chosen to yield a finite limit although this requirement is not forced on us. We will see in the 
following that it is indeed determined to match the three-loop amplitude threshold functions 

7 In the next section we will confirm by an analysis of their Fourier coefficients that they contribute to the two 

different V 6 f? 4 invariants. 
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in ten dimensions. The string perturbation limits of the two series are given by (g = e^ is the 
effective string coupling in D = 6) 

, _ e «4-2e) 

+S w^y [»».-.] 

7-2e ; 


_ —4+e rp 

ll/l | 0 0 0 4- e l _ff f 4 - e0 0 


E 

E 


(4.7a) 


L^v 9 


e E 


oof 


-2 ^( 5 + 2e ) p 

9 f (7 + 2e) [o 3+, ° ] 


, 2e £(2 + 2e)g(4 + 2e)g(7 + 4e) 

9 e(5 + 2e)e(7 + 2e)C(8 + 4e) [oof 


(4.7b) 


The poles at e —>• 0 are contained in the last terms in both lines (corresponding to three loops 
after the Weyl rescaling by g 4 ). The singular terms are explicitly given by 

1 £(3) 


2e £(4)£(6) 


+ 0 ( 6 U ) , 


E 


m 


. 3+€ 


1 


2e ^(4)^(6) 


+ 0(e°) . 


(4.8a) 

(4.8b) 


Note that these two functions are indeed related by T-duality parity. Demanding T-duality 
invariance therefore fixes the relative coefficient between the two Eisenstein series, as anticipated 
in (4.6). From this one can see that the combination given in (4.6) is regular. The non-analytic 
term in g can be deduced from this to be 


nj ( 6 ) 


£(3) 


, , , , lim 

mm 


~9 e + 9 
2e 


2e 


+ 


m 


_ 3 

“ 2 a '^(6)e(8) 


log g T 


(4-9) 


where the adjoint series contributes twice the amount of the spinor series (due to the g 2e 
in (4.7b)). We note that the perturbative limit of the series i7[o|+eo°J also contains a term 
with g~‘ that is formally at loop order L = —1/2. This does not matter for determining the 
logarithmic divergence since this term is cancelled by a corresponding contribution from the 
particular solution as discussed in the introduction to this section. 

The value of the coefficient a can be fixed by comparison with the ten-dimensional three-loop 
correction in the decompactification limit [8, 34, 35]. In a first step, one decompactifies the 
three-loop terms in (4.7) to D = 7 to obtain 



— v 3 € E [o o 3 _ e ] + v 


— v 3+e E [ 3+e o o] + v 


e £(3 — 2e) 

£(6 — 2e) [2-e ° 01 ’ 

-e £(3 + 2e) 
e(6 + 2 e )^ [0 0 2+e] ' 


(4.10a) 

(4.10b) 


Here, v is related to the radius of the decompactifying circle from D = 6 to D = 7 measured 
in seven-dimensional Planck units. Diagrammatically, the expansion above uses the node in the 
lower right corner of the Dynkin diagram, i.e. the vector representation node in the string 
theory convention. The series above are still divergent in the second term but we are interested 
in the first terms that decompactify nicely to give 

%(M) (-^[o o 3] + E [3 o o]) + • • • (4.11) 
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The combinations that decompactifies correctly to D = 10 is in our conventions [19] 


and therefore 


4C(6) 

27 


(E [0 0 3]+ ^[3 0 0]) 


= 406)08) 
2704) 


Putting this together with (4.9) means 


^(o?D^ 503)log<? + ... 


(4.12) 


(4.13) 


(4.14) 


which is indeed consistent with the explicit 3-loop divergence computed in [20], as exhibited in 
[18] through the analysis displayed in [19]. This suffices to determine the value of the anomalous 
term in (1.2c) by acting with the 50(5, 5) Laplacian on the three-loop terms contained in (4.6). 
Acting with the Laplace operator (3.12) on the three-loop terms in 77)0) proceeds through the 
introduction of regularised 50(4,4) Eisenstein series by 


Then 





_i m | a 

2e ((4)((6) ["» S 

1 « 3 > I r 

2e mm [»“*] 


j + O(e), 
+ O(e). 


a d 4 



0 1 
3 J 


3^(3) 

mm' 


(4.15a) 

(4.15b) 

(4.16) 


and the other 50(4,4) series works similarly. Acting with the 50(5,5) Laplacian (3.12) (with 
g = e^) thus yields 


A d 5 T-L 


(6) — A 
(0,1) “ ^ 


£>5 


5C(3) log g + 


4C(6) 

27 


E 


= 20C(3) + 
= 40C(3). 


4C(6) 6g(3) 

27^(4) m 



(4.17) 


This confirms the coefficient in (1.2c). The final function giving rise to the divergent term in 
0 = 6 is then 




( ° ;i) 27£(4) e^o n°°°4- e ] + e(2)^[° 5-*°! 

4C(6)C(8) 


27^(4) 


r ° oo 4 ] + m E [° 5 °o 


(4.18) 
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4.2 Adjoint Eisenstein series of E 8 ( 8 ) and V 6 i? 4 in three dimensions 

Even though the V 6 i? 4 threshold function is regular in three dimensions, it is convenient to start 
our analysis at D = 3, because there is a unique V 6 i? 4 supersymmetry invariant [12]. The exact 
threshold function must satisfy the fourth order differential equation 

{VY i]jk r V){VY lpq]r V)£^ = +6 ib {vr klpq] v) , (4.19) 

consistently with the inhomogeneous equation [15] 

= -198£[ 0 % - {£$ o y . (4.20) 

The tensorial equation (4.19) implies that the wave-front set is associated to the |-BPS nilpotent 
orbit corresponding to extremal black holes in four dimensions (see for example [36]). We 
compute in Appendix A that the automorphic solution to the homogeneous equation 

(VT l[jk r V) {VY lvq]r V)nf^ = 150% (VT klpq p)HZ) , (4-21) 

is proportional to the E 8 ( 8 ) Eisenstein series 


Ar o 
oooooo 



(4.22) 


This is in some sense the ancestor of all three-loop terms in higher dimensions. For this ‘adjoint’ 
A 8 ( 8 ) series the Laplace eigenvalue is 


AEr 0 1 = 2s(2s - 29)Ar 0 i. (4.23) 

[o 0 0 0 0 0 s] [o 0 0 0 0 0 s] 

The function is regular at s = -y. 

We can perform the decompactification of this function from D = 3 to D = 6 by computing 
the constant term in the parabolic subgroup associated to the 6 th node of A 8 ( 8 ) with Levi factor 
Lq = GL( 1) x Spin(5,5) x SL(3). For arbitrary s the result is 



= 3 ^(^ ,s ~ H)£(2s — 12)g(2s — 13) ,425+11/4 ^-3) 

£(2s)£(2s - 5)£(2s - 9) [o^«S] £(2s) [ 000 s -°f 

q 5 _ a)/2 Z(2s - m(2s - 19)g(2s - 20)g(4s - 29) 

£(2s)£(2s - 5)£(2s - 9)£(4s - 28) [° ° o s 0 9 ] + ''' ’ 


(4.24) 


where we have only listed the SL(3)-singlet contributions. The GL( 1) parameter v here is related 
to the volume of the decompactifying three-torus. The last term in (4.24) is subdominant for 
s —» ^ and the first term can be rewritten according to the functional relation 


_ £(2s - 4)£(2s - 5)g(2s - 6)g(4s -7) 

[° • 0 S ] £(2s)£(2,s — l)£(2s — 2)£(4s — 6) [° i- o ° ] ' 


(4.25) 


Then the dominant pieces become 


3 g(2s - 15)g(2s - 16)g(2s - 17)g(4s - 29) ^ 

V £(2s)£(2s - 5)£(2s - 9)f (4 s - 28) 


| (28+l)/i ^( 2s 3 ) E 

«2 8 ) [» 



(4.26) 
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that tend to 


v 


3 


m (t(5)p 

e(ll)U(2) [°i 


+ E 


(4.27) 


for s —^ - 7 j-. This is precisely the combination appearing in the ansatz (4.6) in six dimensions. 
From this we conclude that the D = 3 ancestor of the correct D = 6 three-loop divergence is 
given by 


0/(3) 

'*'( 0 , 1 ) 


“ «s) E 


0 0 0 0 0 0 ^ 


4C(6)g(ll) 

27£(4) [oo o ooo y] 


(4.28) 


Here, we have defined the normalisation constant b. We will now use this finite adjoint Eg 
function to determine the remaining terms in (1.2c) and (1.3). 


4.3 Form factor divergence in D = 5 

We consider the decompactification limit of the general adjoint E 8 ^ function from D = 3 to 
D = 5. The Levi subgroup in this case is L 7 = GL{ 1) x SL{ 2) x E§t§y One has 

_ . 4 e(2a - 12)e(2a - 11) (2s+1)/3 g( 2 s - 2 ) 

£(2s)£(2s-5) [ 00 S 7oo] e(2s) [ 00 S 0 .- 1 ] 

(4.29) 

where we have focussed on the pieces that are relevant for the discussion. These are in particular 
singlets under the SL{ 2 ) group of the decompactifying two-torus whose size is related to the 
parameter v. The first term in the above expression is clearly divergent for s —> ^ since £(0) 
diverges and the adjoint Eq^ function is regular at s = 1. We leave the divergence implicit in 
the second term in (4.29). 

We will be interested in the string perturbation limit of the individual terms. For the adjoint 
E$r 6 ) series one finds from (2.13) 
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(4.30) 


where these are now functions on the symmetric space SO(5,5)/(SO(5) x SO(5)) parametrised 
by the NS moduli in five dimensions, with the standard H 5 labelling and g is the five-dimensional 
string coupling. According to (4.29) we need this expression in the limit s —> 1 + e. In this limit 
the last term disappears and we are left with (recalling the £ prefactors from (4.29)) 


g(2)£(2e) 

£(6)£(11) 



+ 9 


-4 £(3) 
£( 2 ) 


+ 0(e) 


(4.31) 
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The string perturbation expansion of the fundamental -E?6(6) series in (4.29) is 
_ rr —8(6 —^~)/3 C(2g — ll)£(2s — 8) ^ g-(l'o- 2 s)/ 2 , i(?‘ s ~ 5) 
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= 9 
+ g s E 

which exhibits the expected pole when s —>• | 

More precisely, for 

C(9)/C(H) from (4.29)) 


£(2s) 0 ° a -f 


(4.32) 


(which is the right value after taking into account 

9 
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the shift from above). More precisely, for s = | + e we have the terms (recalling the finite 
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(4.33) 

(4.34) 


£(6)£(n) e(6)e(n) 

where we have used a functional equation in the second step. 

We collect all the relevant terms from above and find for s = + e (ignoring the u 4 that 

only indicates the decompactification) 
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(4.35) 


Here, we have used that £[o o o ° ] = £|oooJJ (for this particular value of s = 1) as required by 
T-duality. Multiplying by the normalisation b from (4.28) one has then as for the logarithmic 
term in g: 
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where we have used that the £g(6) Eisenstein series appearing in the R 4 corrections has the 
string perturbation expansion 

-2<£(2) 


/ Er o 1 

ful Ll 0 0 o oj 


= 5 3 + <L 


e(3) 


E 


(4.37) 


and £(om) = 2C(3)£[a 0 ° 00 ]. ^ 4e coefficient in (4.36) matches the claimed coefficient in (1.3). 

The next task is to determine the coefficient in the Laplace equation. This one gets by 
acting with the Laplacian on the regularised series. To this end we note the following functional 
relations 

_ f (2 a - 19)f(2s - 17)£(2s - 16)f(4s - 29) 
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This absorbs the pole in the first term in (4.29), leading to 
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Both Eisenstein series are singular and need to be regularised. The Laplace eigenvalues for 
s = T e are 


A E 

A E 


2“ e 

0 0 0 0 0 


0 0 0 0 S+e 


= 2(1 + e)(2e-9)£Y a 


2 e 
0 0 0 0 0 


= ( -e" + 8e - 18 ] E 


(4.40) 

(4.41) 


0 0 0 0 2+e I 


and from the analysis of the poles above we know that we can define the regularised series 
(denoted with a hat) by 
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(4.43) 


Therefore we find the anomalous term in the Laplace equation for £m,i) t° be 
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(4.44) 


Here, we have used the normalisation b from (4.28). This matches the claimed coefficient 
in (1.2c). 


4.4 Form factor divergence in D = 4 


We need to consider the decompactification of the adjoint E S ( S \ series to D = 4. The Levi 
subgroup is now Lg = G'L(l) x E 7 (jy The constant term formula (2.13) leads to a number of 
terms of which we only display the ones relevant for the derivation of the form factor divergence: 
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Both coefficients are regular but the two Eisenstein series diverge at s = 4j-. The parameter r 
here is related to size of the decompactifying circle. 

Next, we require the string perturbation limit of the two Eisenstein series on Ej^ / (SU (8)/Z 2 ). 
For the adjoint E 7 ( 7 ) series one has 
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We have chosen a representative that brings out the divergence at s = if + e explicitly through 
the factor £(2s — 11) = £(2e) = — ^ For s = ^ + e one then has (after reinstating the 

prefactor from (4.45) and not exhibiting the r dependence) 
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where we have used the string perturbative expansion 
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For the fundamental ^ 7 ( 7 ) series one finds the perturbative string expansion 
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This expansion is exact and we have again chosen a form that brings out the divergence explicitly 
through the prefactor £(2s — 10) = £(2e) = + ... for s = ^ + e. Reinstating the prefactor 

from (4.45) one then obtains for the j and log g terms 
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Putting (4.47) and (4.50) together we see that the divergence indeed cancels out and the 
logs terms are 
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The reassembly of the full adjoint 74 ( 7 ) series was only done at the level of the string perturbative 
terms but will hold fully. Multiplying in the normalisation b from (4.28) the final result is 8 
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s Recall that f(i 0 ) = C(5)F^5 0 q 0 0 0 ] without a factor of 2. 
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This confirms the relevant term in (1.3). 

We now turn to the computation of the anomalous term in the Laplace equation. We define 
regularisations (denoted by hats) of the series appearing in (4.45) in accordance with the pole 
analysis at s = + e by 
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+ 0(e). (4.54) 


Here we have kept the prefactors as they appear in the combination (4.45). Using now 
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The anomalous term in the Laplace equation of } is therefore 


(A + 60) — b ^7 + 

thus confirming (1.2c). 
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5 Fourier expansions and Whittaker vectors for 50(5, 5) 

The Fourier coefficients of the functions £{p^ q) are constrained by supersymmetry [16, 17, 10]. 
This can be rephrased in terms of nilpotent orbits [17] and constraints on degenerate Whittaker 
vectors [37, 38]. In this section we will investigate these issues for the case of D = 6 where the 
Cremmer-Julia group is 50(5,5). 


5.1 Fourier modes 

Given a U-duality invariant function f(Vo) on G/K, as for example one defines the 

(abelian) Fourier coefficients in some (abelian) unipotent subgroup U C G by 

Fmu(VD) = [ /{uVD^muiujdu. (5.1) 

Ju 1 

Here we have used the following notation. 
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• U 1 = U (Z)\[7 where U (Z) = U HG(Z) denotes the U-duality shifts contained in the chosen 
unipotent U. In an appropriate normalisation the integral is nothing but an integral over 

[ 0 , 



• The unipotent subgroup U corresponds to a choice of subset of moduli in D dimension. All 
these moduli must be axionic, i.e. , have shift symmetries, and their shift symmetries must 
close on U. If U is abelian then all shifts commute but more general cases are possible, 
for example in the presence of NS5-branes [39, 40]. 


• As for abelian U all shifts commute we can diagonalise them simultaneously and denote the 
corresponding eigenvalues by mjj £ Z dim ( u K Similarly, we can parametrise any element 
u £ U by u = exp (xu’tu) where tjj denotes a basis of appropriately normalised generators 
of the unipotent group U. The character (u) appearing in the formula above then is 
nothing but 

VW H = exp (2 ttm% ■ xu) , (5.2) 

i.e. a collection of phases associated with the axionic moduli xu that are physically defined 
modulo an integral shift. 


• If U is non-abelian, then only the abelianised part is captured by the integral (5.1). When 
U is non-abelian the derived series of U defined by U ^ = \U^ k ~ l \ U^ k ~ 1 ^] (with C/(°) = U) 
is non-trivial for 0 < k < n for some n > 0. One can then define similar non-abelian Fourier 
coefficients for each of the U ^ for k > 0. 

The integers my correspond physically to instanton charges and the Fourier coefficients arrange 
themselves on orbits of the action on U of the (reductive) Levi subgroup L associated with the 
unipotent U. The case rhy = 0 corresponds to i= 1 and therefore the zero mode Fourier 
integral reduces to the constant term integral (2.13) discussed in section 2. 


5.2 Fourier modes of SO(5,5) series 


We consider the Fourier modes in an expansion corresponding to decompactihcation to type IIB. 
At the level of Lie algebras this means that we are looking at the decomposition of so (5, 5) under 
Lie(Z/ 2 ) = 0l(l) ©sl(2) ©s((4) where sl(4) = so(3,3) describes the decompactifying four-torus, 
fll(l) is related to its volume and s((2) is the type IIB S-duality. The relevant decomposition is 

so(5, 5) “ 1 ( “ 2) © (2, 6) (-1) © (fll(l) © sl(2) © s[(4)) (0) © (2, 6) (1) © 1 (2) (5.3) 

'-—V-' 

Lie (Us) 


The superscripts in this equation denote the 0((1) weight and the representations of sl(2) ©sl(4) 
are labelled by their dimension. As indicated, the Lie algebra of the non-abelian unipotent 
subgroup U 2 in this expansion consists of two pieces occurring at weights +1 and +2 with respect 
to 0l(l). The associated Heisenberg algebra is realised on the symmetric space SO( 5,5)/(50(5) x 
50(5)) through the Killing vectors 
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(5.4) 




satisfying 


(5.5) 


K,4 l ] = \s a ^ kl h. 

The indices i,j = 1,... ,4 here are fundamental sl(4) indices and an antisymmetric pair [ij] cor¬ 
responds to the vector of so(3, 3) = s((4). The index a = 1, 2 is a fundamental sl(2) index. An 
abelian Fourier mode defined in (5.1) here is labelled by an SL( 2,Z) doublet Q = (p,q) = mu 
of SO( 3,3) vectors. The non-abelian Fourier mode associated with the weight +2 singlet corre¬ 
sponds to a single integer N and we will write the most general Fourier mode as FQ^iVe). Note 
that this notation is not uniquely defined because not all Kq can be diagonalised simultaneously 
when N ^ 0. A related discussion can be found in [40]. For a purely abelian Fourier mode Fq q, 
the action of the generators is ^Fq^Vq) = 0 and k^Fq^q (Vg) = 2niQa FQfi(V&). A non-abelian 
Fourier coefficient Fqjg satisfies in general k^FQ^fVe) = 2niN FQ^(y&)- In physical terms, the 
Fourier modes Fq jv correspond to fundamental and D1 strings wrapping the T 4 and define Qa 
whereas the additional single integer N gives the number of Euclidean D3 branes wrapping T 4 . 

The space of vector doublet charges (p, q ) € Z 12 stratifies under the action of 5L(2, K.) x 
50(3,3) according to the number of linearly independent vectors spanned by p and q and the 
sign of their norm. This stratification is isomorphic to the stratification of real nilpotent orbits 
of 50(5,5) of dimension smaller than 26, and is displayed in figure 1, where the number of 
+, —, 0 gives the number of linearly independent vectors whose norm is respectively positive, 
negative or null. 9 * * 
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\ / \ / 

h,0] [+,0] 
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{ 0 } 


(o,i: 
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1 1_' 
^2’ 2> 

I 

( 1 , 1 ) 


(3.0) 


Figure 1: Closure diagram of the real nilpotent orbits of 50(5, 5) of dimension smaller than 26. The 
second copy exhibits the fraction of supersymmetry charges of the two chiralities preserved by the corre¬ 
sponding instanton. The total fraction of supercharges that are preserved is given by half the sum of the 
two chiral pieces. 

9 In terms of signed partitions of 10 parametrising the nilpotent orbits of 50(5, 5), a + and a — correspond 

respectively to a 3-box line with respectively -|-h and —I—, whereas a 0 corresponds to a doublet of 2-box 

lines, the remaining being understood as being 1-box lines with the appropriate signs to neutralise the partition, 

see e.g. [36]. 
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Condition on (p, q) £ Z 12 

‘Dimension’ of 

part of nilpotent 50(10, C) orbit O 


subspace 

Bala-Carter label 

weighted Dynkin 

dime O 



of O 

diagram of O 


bl 2 lb 2 -(p'?) 2 /o 

12 

A 2 (special) 

[° 2 0 0 ] 

26 

\p\ 2 \q\ 2 - (p ■ q) 2 = 0 

11 

2>A\ (not special) 


24 

bl 2 = bl 2 = (p- q ) 2 = 0 

9 

{2Ai)" (special) 

f° 0 0 1 

20 

p A q = 0 

7 

{2Ai)' (special) 

to 

0 

0 

0 0 

16 

b| 2 = \q\ 2 = p A q = 0 

6 

Ai (special) 

0 0 

0 

0 

14 


Table 1: Orbits of 50(3, 3, Z) x SL( 2, Z) acting on the twelve charges (p,q) £ Z 12 describing the 
(abelian) Fourier coefficients associated to P 2 . The intersection with the complex nilpotent orbits of 
50(5,5, C) = 50(10, C) is also given. 


The algebraic constraints defining this stratification are displayed in table 1, where we have 
also indicated the nilpotent 50(10, C) orbits that are intersected by the orbits of the 50(3,3) x 
SL(2) action on 1R 12 and have labelled them both by their Bala-Carter type and their weighted 
Dynkin diagram [41, 42, 43]. The qualifier ‘special’ or ‘not special’ refers to Lusztig’s property of 
nilpotent orbits [31] and is related to whether they can arise as wave-front sets of automorphic 
representations [32, 24]. The information of this table can partially be extracted from [44, Sec. 
5.5.2]. The notation here is such that \p\ 2 = p ■ p denotes the 50(3, 3)-norm of the vector p and 
(p ■ q) 2 denotes the square of the 50(3,3) inner product. The wedge p A q is a tensorial object, 
namely the outer product of the two vectors. 

The non-trivial (abelian) Fourier coefficients must fall into the topological closure of one 
of the above classes. They correspond to space-time instantons that are defined in the back¬ 
ground of supersymmetric solutions in supergravity. These solutions are defined in the Euclidean 
signature as solutions to the pseudo-Riemannian non-linear sigma model over G/K* (for a non¬ 
compact real form K* of K) satisfying the energy momentum tensor constraint 

{P^),P v {^))d^ ®d^ v = 0 . (5.6) 

Such a spherically symmetric solution can be obtained for any representative p of the nilpotent 
orbit restricted to the coset component g © £*, as [45] 

V($) = exp(—, (5.7) 

and the corresponding solution preserves a given amount of supercharges depending of the com¬ 
plex K( C) orbit of p, as is displayed in figure 1. The smaller the orbit is the stronger the 
constraints from supersymmetry are, and figure 1 encompasses all the BPS orbits. The differen¬ 
tial constraints satisfied by the BPS protected threshold functions are themselves associated to 
harmonic superspace constructions of the associated linearised supersymmetry invariant as su¬ 
perspace integrals of BPS protected integrands preserving the same amount of supersymmetry, 
i.e. 1/2 BPS for i? 4 , 1/4 for V 4 i? 4 and 1/8 BPS for V 6 i? 4 [10], 
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For a generic adjoint Eisenstein series 


E 


(5.8) 


the Fourier coefficients are restricted to the BPS orbits displayed in figure 1, and generically 
cover all of them non-trivially, but for special values of s for which they are further restricted. 
The constraints on the Fourier modes follow from the tensorial differential equations the corre¬ 
sponding function satisfy [10, 12]. One can check explicitly that these differential equations are 
satisfied by the generating character of the series, but one must be careful when the Eisenstein 
series is not absolutely convergent. Indeed we will see in the next subsection, that although the 
adjoint Eisenstein series at s = | is generated by a character satisfying the quartic constraint 
associated to the dimension 24 nilpotent orbit of type 2>A\ in table 1, its Fourier modes do not 
satisfy the corresponding constraint |p| 2 |qj 2 — ip,q) 2 = 0. This nilpotent orbit is indeed not 
special, and so one understands this property from the analysis of [32, 24], that states that the 
maximal orbits in a wave-front set of an automorphic function are special. There are only two 
non-special nilpotent orbits of 50(5,5, C), as we display in figure 2. 

From the tensorial differential equation one derives that 
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\p\ 2 \q\ 2 ~(p,q) 2 ^0 (m)£2 12 


bl 2 = l?l 2 = (p,g) = o 

p A q = 0 

M 2 =pAq= | q\ 2 = 0 


(P, q) e z 12 n O 9 
{p, q) G z 12 n o 7 
(p, q ) G z 12 n o 6 


(5.9) 


where O" is the dimension n algebraic variety in R 12 of solutions to the corresponding constraint. 
Here the value of s is understood to be generic, whereas the support of the Fourier modes reduces 
at some specific values 
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(5.10) 


where the differential equations degenerate. Strictly speaking the undisplayed coefficients vanish 
as an artefact of Langlands normalisation, but we will see in the next subsection that these 
relations are indeed satisfied in a suitable normalisation in which the coefficients are finite. The 
value of N is never constrained. When N = 0, one can define the abelian Fourier modes ip, q). 
and the reduction of a differential operator in a given representation implies the same constraint 
as the corresponding algebraic equation. We will moreover assume that the dimension of the 
algebraic variety that supports the non-abelian Fourier modes is still correctly determined by 
the algebraic solution to the constraint to be half the dimension of the corresponding nilpotent 
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Figure 2: Closure diagram of the nilpotent orbits of S'O(10, C) exhibiting the special orbits. The labels 
on the nodes represent partitions of 10 and have to be completed by Is as needed, such that for example 
the label 3 corresponds to the partition 3-1'. The two orbits in small font, 3 • 2 2 • l 3 and 5 • 2 2 • 1, are 
not special. Removing them together with lines connected to them one obtains an ‘up-down’ symmetric 
diagram as required by the Spaltenstein map. 


orbit, even if the structure of the solution is much more complicated in general for non-abelian 
Fourier modes. 

Let us explain this in more detail in some examples. The vector Eisenstein series satisfies in 
general 


Dfg Er o ' 
iD Loo” 


s(s — 4) 


li6 E 


(5.11) 


and using the 3-graded decompositions associated to (5.3) 


16 ^ 4 (-1) © ( 2 , 4 ) (0) © 4 (1) , 

10 ^ 2 ( ~ 1) © 6 (0) © 2 (1) , (5.12) 


one computes that the restriction of (5.11) to the degree 2 component gives 

-£ ipk l£ afi Q^Q 1 fF Qfl ^ 0 0 o j = 0 => q Ap = 0 . 


(5.13) 
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Similarly for the spinor Eisenstein series 
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one obtains the degree 2 component 
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(5.14) 


(5.15) 


More generally Eisenstein series satisfy tensorial differential equations that are associated to the 
algebraic constraints defining the corresponding nilpotent orbits. Because the maximal degree 
component of these differential equations in a given parabolic gauge only involves the axion 
fields, they become algebraic constraints in Fourier space that reproduce the corresponding 
constraint. One can in this way understand the wave-front set from the differential equations 
satisfied by the Eisenstein series [12]. Strictly speaking the reduction of (5.14) to the algebraic 
constraint (5.15) on the Fourier modes (p,q) is only valid when N = 0, whereas in general one 
gets a corresponding differential equation for the axion with a right-hand-side linear in N that 
depends explicitly on the axion moduli. Nonetheless, the constraint still makes sense for counting 
the lattice dimension, i.e. the dimension n of the algebraic subvariety O" in the corresponding 
vector space of charges IR 12 such that the Fourier modes are valued in Z 12 D O". We can see this 
explicitly in the perturbative string theory limit, in which one decomposes instead the series in 
RR Fourier modes e 2m ( p ’ a \ with the Spin(4,4) Weyl spinor axion a. We compute in Appendix 
B using the Poisson summation formula that 
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0 0 0 
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-sp .-4«2 (s-2)) 
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(p,p)=0 
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n|gcd(p) 


2(2—s) A 9 


\Z{p)\ 


*V10 0](^)A^_2(27r 5 - 1 |Z(p)|) e 2 ^’ a ) (5.16) 


where v p is the SL{4 ) subgroup of 50(4,4) that stabilises the null vector p , and \Z(p)\ 2 is 
the 50(4,4) moduli dependent invariant mass associated to the charge p. In this equation 
one finds that the RR Fourier modes are constrained to be null vectors, and the three extra 
integer sums come from the 5L(4) Fourier modes of the Eisenstein series in the fundamental 
representation, with total ‘lattice dimension’ 10, as for the ten-dimensional lattice associated to 
the spinor Eisenstein series in the M-theory limit. In this case however, we have a non-abelian 
decomposition for which each null vector p of Z 8 defines a particular Fourier decomposition of 
the corresponding SL( 4) C 50(4,4) moduli space. The graded decomposition of the vector 
representation indeed implies that for a generic charge 
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one has the quadratic constraint 


Qio = 


/ 0 
0 

V 0 


Qs p 

Qi 
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(p,p) ^ 
Qs p 

0 / 


= 0 , 


(5.18) 


such that the £[ s .i 00 ](wp) must satisfy the quadratic differential equation associated to the 
minimal nilpotent orbit 


Dg £[_i o 0] 


(s - l)(s - 3) 


l6^[ s -i o o] 


(5.19) 


For s = 1, the ST(4) Eisenstein series is a constant, and one gets back T|§ o o ° J. 

Let us now consider the inhomogeneous equations (3.2,3.3). The Fourier modes of £ (0i o) of 
vanishing D3-brane charge are defined by doublets of proportional null vectors and so its square 
involves the generic sum of doublets ( pi , qi). One computes that the sum of such doublets gives 
generically a doublet of linearly independent non-null vectors (p, q) of opposite signature, that 
has a negative quartic SL( 2) x 50(3,3) invariant 


bl 2 M 2 - (p,q) 2 < o 


(5.20) 


This can easily be seen if we take 


(pi,Qi) = (p,0) , (P 2 ,q 2 ) = (0,q) , {p,q)^ 0. (5.21) 

This corresponds precisely to the structure of the Fourier modes of a generic function E\o a o ° j, 
although its Fourier modes can have the two signs for |p| 2 |g| 2 — (p, q ) 2 in general. To see this one 
uses the property that the ‘fundamental’ representations are all 3-graded in this decomposition 
(5.12) such that the positive degree generator in (5.3) to the third power vanishes automatically 
in all these representations without any further restriction. Note that the real 50(5,5) nilpotent 
orbit associated to instanton charges of negative quartic invariant corresponds to a complex 
Sp( 4, C) x Sp( 4, C) orbit of a non-BPS solution through the Kostant-Sekiguchi correspondence 
[43] . Nonetheless, a given representative does not need to lie in the intersection of the respective 
orbits that are related by the Kostant-Sekiguchi isomorphism. This additional restriction only 
applies when one requires moreover the corresponding black hole solution to be regular in seven 
dimensions [36], but should not apply to the space-time instanton as such. 

One can also check the representation of the Fourier modes in the parabolic subgroup asso¬ 
ciated to the decompactification to eight dimensions 

so(5,5) 9* 3 (-4) ®(2®2® 3) (-2) 0 ( fl t x 0 sl 2 0 sl 2 © s( 3 ) (0) © (2 ® 2 ® 3) (2) 0 3 (4) , (5.22) 

with 


16 9 * 2 ( ~ 3) 0 (2 ® 3 ) (-1) 0 (2 ® 3 ) (1) © 2 (3) , 
10 9 * 3 ( ~ 2) 0 ( 2 ® 2 ) (0) 0 3 (2) . 
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The Fourier modes give a triplet of 50(2,2) vectors pi associated to the effective Euclidean 
1-brane coupled to the 6 vector fields along T 2 , and a conjugate triplet of singlets q l associated 
to the effective Euclidean 2-brane coupled to the 3 tensor fields. 10 For the function E\o s o ° j, the 
condition that the differential operator in the spinor representation reduces to the third order 
implies 

s ljk Pi A pj A p k = 0 , (5.24) 

such that only two of the vectors are linearly independent. The counting of the number of 
modes here works as follows. There are two vectors (4 + 4) plus a third depending on the two 
others (+2) and the three weight 2 components (+3), and we get back a 13 dimensional lattice 
of Fourier modes. 

For the function E^o o o ° J, we moreover require the vector representation differential operator 
to reduce at second order, such that all the scalar products vanish 

(Pi,Pj) = 0, (5.25) 

getting therefore three more constraints (since ps is already assumed to be a linear combination 
of pi and P 2 ), recovering the 10 dimensional lattice of Fourier modes. 

The minimal representation corresponds to the restriction in which all pi are proportional 
and null, and moreover piq 1 = 0, such that one gets a3+l+l+2=7 dimensional lattice 
of Fourier modes as required. The generic sum of two such vectors gives a doublet of linearly 
independent vectors, but cannot give three linearly independent vectors, such that one gets 
back that the source term (3.2,3.3) indeed sources generic Fourier modes of a function of type 
F^oso°|, and is consistent with it in the sense that it does not source more generic Fourier 
modes. 


5.3 Whittaker vectors of 50(5,5) series 

In this section we analyse the degenerate Whittaker vectors of some maximal parabolic 50(5, 5) 
Eisenstein series. Whittaker vectors are special cases of Fourier coefficients when the unipotent is 
taken to be maximal unipotent N of all unipotent elements in a Borel subgroup. There is a close 
connection between Whittaker vectors and Fourier modes associated to nilpotent orbits [23, 38]. 

The number of instanton charges that label Whittaker vectors is equal to the real rank of 
the group. For 50(5,5) we therefore have five instanton charges rfiN that can be arranged on 
the Dynkin diagram of 50(5, 5). When some of the charges vanish, Whittaker vectors are called 
degenerate and a general formalism for determining degenerate Whittaker vectors was presented 
in [37]. We will write a Whittaker vector as 


W 1 
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m 4 I _ 


N(Z)\N(R) 


Ex „4l^r m 4 I j 

I "1 s 3 S5 J [™1 m 2 m 3 ^ I 


(5.26) 


10 By effective Euclidean p-brane we mean any Euclidean p+fc-brane that wrap p directions of the decompactified 
T 2 and k directions in the other T 2 . 
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with iJj ^ m\ m -2 m 3 "‘j* | defined as in (5.2). The Whittaker function types can be labelled by the same 
labels as the (complex) nilpotent orbits and we will refer to them by these names. The wave- 
front set of an Eisenstein series will be the largest complex orbit with non-vanishing Whittaker 
vectors, which is unique [46, 32]. 11 In the following we will always evaluate the Whittaker vectors 
at the origin of moduli space, i.e. the identity of the Cartan subgroup. 


5.3.1 50(5,5) vector series 


The series 


g(2s-6)g(2s — 4) 
£(2s)£(2s-3) 


(5.27) 


can have Fourier coefficients at most associated with the orbit of type {2Ai)' since it can be 
functionally realised on the coset Pi\50(5,5) of dimension 8 . The functional dimension of an 
automorphic realisation is half the dimension of the maximal nilpotent orbit contributing to 
the Fourier coefficients and therefore inspection of table 1 shows that there can be no Fourier 
coefficients beyond the orbit (2A\)'. 

The Whittaker vectors corresponding to the {2A\)' orbit can be represented by instanton 
charges with Dynkin diagram |^ooo™J for m,n G Z. By contrast, the instanton charges for the 
(2A\)" orbit are for example of the type 0 n ° J. One can check that the associated Whittaker 
vectors vanish for all values of s in (5.27). For a generic s in (5.27), the Whittaker vector can 
be evaluated with the help of [37] to be 


W, 


s 0 0 


0 1 
0 J 



g(2* ~ 3) 

£( 2 s) 


W [s-z] (M) W [s -§] (M) > 


(5.28) 


where we have placed ourselves at the origin in moduli space for simplicity. Here, the function 


W [s] ([n}) 



2s )7^ s _i/ 2 (2vr|n|) 

k\n 


(5.29) 


is the Whittaker function for the SL( 2) series E\ s 1 at the origin of moduli space. It vanishes 
linearly for s —> 0 (and s —> |) due to the £(2 s) denominator. Inspecting (5.28) then shows that 
the whole Whittaker vector for the 50(5,5) vector series vanishes for s —> 

Performing a similar calculation for the smaller A\ orbit shows that the associated Whittaker 
vector never vanishes (as a function of s 7 ^ 0) confirming the reduction of the Fourier coefficients 
already mentioned above. 


5.3.2 50(5,5) adjoint series 

The adjoint Eisenstein series of 50(5,5) 


E 


f (2s - 6 )f (2 s - 5)£(2s - 4)^(4s - 7) 
e(2sK(2s-2)e(2 S -l)^(4s-6) 


E 


(5.30) 


11 i.e. the closure of this orbit includes all the orbits for which the Whittaker vectors do not vanish. 
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can be realised as a function on the space P2\SO(5,5 ) of dimension 13. The maximal nilpotent 
orbit supported by this function for generic s is therefore the orbit of dimension 26, labelled 
type A 2 in table 1. 

For special values of s one has reductions of the orbit type. As indicated in (5.30), the 
adjoint series has a functional relation that relates s -H- ^ — s. In Langlands normalisation, one 
can use the Langlands constant term formula to show that E^o s 0 ° J has simple zeroes at s = 
s = 1 and s = §. Inspecting the prefactor in the functional relation one concludes that there 
has to be a simple pole for s = 2, s = | and s = 3. These have to be taken into account when 
discussing the simplifications in the Whittaker vectors. 

Performing this analysis here implies that the degenerate Whittaker vector of type A 2 takes 
the value 


W 1 


g(2s-3)£(2s-4) 

£(2s)£(2s - 2) 


W[s-l/2 s—2] 


([HI), 


(5.31) 


for unit charges. Let us discuss various limits of this formula. 


• In the limit s —>• the 5L(3) = A 2 Eisenstein series becomes of minimal type [40] and the 

corresponding Whittaker vector vanishes linearly as does the prefactor. This means that 
the suitably normalised adjoint 50(5, 5) Eisenstein series has no A 2 Whittaker vectors for 
s = This is consistent with the fact that it is related to a different Eisenstein series 
through 
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0 s 0 


0 1 
0 J 


£(2s -2) 

f(2 s) [l-o-* 


(5.32) 


This shows that in the limit s —y ^ the suitably normalised adjoint series is proportional to 
the vector series at s = 1 and therefore has Whittaker vectors of type (2Ai)' (and smaller). 


• In the limit s —> 1, the prefactor in (5.31) vanishes linearly and so does the Whittaker 
vector on the right-hand side, giving a total vanishing up to quadratic order. In view of the 
linear vanishing of the Eisenstein series in this limit this means that the adjoint 50(5,5) 
series has no A 2 Whittaker vectors in the limit s —>■ 1. Turning to the 3A± Whittaker 
vectors shows that they also vanish in this limit. Further analysis shows that there is an 
effective reduction to type (2A\)" for the Fourier coefficients. The adjoint series indeed 
satisfies the functional relation 


£(2s - 4)£(2s - 3)g(4s - 7) 
f(2s)£(2s - l)f(4s - 6) 


(5.33) 


showing that at s = 1 there is a relation to the chiral spinor series that will be discussed 
below. 


• In the limit s —> |, the SL( 3) = A 2 Eisenstein series becomes of minimal type [40] and the 
corresponding Whittaker vector vanishes linearly as above. But in this case the prefactor 
remains finite. As the whole Eisenstein series vanishes for s —> |, there is in fact no 
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simplification in the degenerate Whittaker vector if one removes the overall vanishing by 
suitable normalising £ factor. If one performs a similar analysis for the 2>A\ type Whittaker 
vectors one also concludes that there is no simplification in this limit. Therefore, the 
suitably normalised s = | adjoint Eisenstein series has Whittaker vectors of all types up 
to A 2 and indeed is of the same type as the generic adjoint Eisenstein series. 

However, the s = § adjoint function differs from the generic case in that the Eisenstein 
series £|o§ 0 °] is square integrable according to Langlands’ criterion [25, §5]. This is 
similar to what happens for the minimal series (for D < 6) and the next-to-minimal 
series (for D < 4) [17] and signals that it belongs to a small unitary representation. The 
generating character of the adjoint Eisenstein series at s = | does satisfy to an additional 
quartic differential constraint associated to the 2>A\ nilpotent orbit, and by Langlands 
formula all its non-vanishing constant terms (in e <^Wp,^(V))) ^ ag we ]p although its 
Fourier modes violate this constraint as we have just exhibited. Moreover, the Laplace 
eigenvalue is —12 and thus outside the range of the associated continuous (degenerate) 
principal series of solutions to (3.31) with s = \ +ir (r € R) that has Laplace eigenvalues 
given by —^ — 4 r 2 . It therefore is part of the discrete spectrum of the Laplace operator. 12 

Note that the Whittaker vector only spans the Fourier modes associated to the non-BPS 
real nilpotent orbit of 50(5,5). The notation A -2 means that a representative of the nilpotent 
orbit can be defined as a linear combination of the simple roots of a subalgebra 5(3 C so(5,5). 
In this case one can indeed realise a representative of the nilpotent orbit of dimension 26 as a 
linear combination of the two simple roots of 3(3 through the embedding 

SL{ 3) C 50(1,1) x 5L(3) x 50(2,2) c 50(3,3) x 50(2,2) c 50(5,5) , (5.34) 

so that the Levi subgroup of the real stabilizer of the corresponding nilpotent orbit is 50(1,1) x 
50(2,2). One straightforwardly works out that this is also the stabilizer of a doublet of linearly 
independent non-null vectors of opposite signature in 5L(2) x 50(3,3). However, a doublet of 
linearly independent non-null vectors of the same signature is stabilized by 50(2) x 50(3,1), 
which is the centralizer subgroup of 517(2,1) in 50(5,5), through the embedding 

51/(2,1) c 50(2) x 51/(2,1) x 50(1,3) C 50(4,2) x 50(1,3) C 50(5,5) . (5.35) 

Nonetheless, it seems clear that the adjoint Eisenstein series should have non-zero modes in all 
these three real orbits, as expected from the discussion in [47]. 


5.3.3 50(5,5) spinor series 

The chiral spinor series 
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000 
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£(2s - 5)£(2s - 7) 
£(2s)£(2s — 2) [000 4-. 
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(5.36) 


12 An analogous phenomenon happens for the adjoint 1 / 7 ( 7 ) series at s = 4 that was mentioned at the end of 
section 2.3. 
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is related to the anti-chiral spinor series as shown. It can be realised through functions on the 
space I 3 4\S'0(5,5) of dimension 10. Referring back to table 1 we conclude that at most the type 
(2A\)" orbit can appear in the Fourier coefficients. The series E|o o o “ J has simple zeroes for 
s = 7} and s = |. 

The type {2A\)" Whittaker vector is 


Wr 


[ oo °; 


£(2 s - 3)£(2 s - 3) 

^(2s)e(2s - 2) 


w [ s -%] (N) w [s-i) (h) 


(5.37) 


Note that this formula is perfectly consistent with (5.16), identifying the Whittaker vector of 
the SL( 4) function 

W[ 3 _i o o] ([o,n,o]) = ^ _ 2 j lF| 3 _3] (W) , (5.38) 

one is left with a remaining factor of Wj- a _3 j ([gcd(p)]) for the integral null vector p that 

matches precisely (5.16). 

We again discuss various limiting values of (5.37) for s. 


• In the limit s —> i, the prefactor in (5.37) vanishes linearly but the two SL( 2) Whit¬ 
taker functions remain finite. This is in agreement with the linear vanishing of the whole 
Eisenstein series and means that a suitably normalised version does not exhibit any sim¬ 
plifications in this limit. 


• In the limit s —> 1, the prefactor vanishes linearly and the two Whittaker functions in (5.37) 
remain finite, implying an overall linear vanishing of the {2A\)" Whittaker function (5.37) 
for the .SO^S) chiral spinor series (since the Eisenstein series itself is regular at this 
value). Indeed, there is a functional relation 


E 


£(2s — 4) 
f(2 s) [§-°° 



(5.39) 


showing that in the limit s —y 1, the chiral spinor series will have the same behaviour as 
the vector series (5.27) at s = |. We already showed in section 5.3.1 above that in this 
case one has a reduction to the minimal A± orbit. 


• In the limit s —>• |, the prefactor in (5.37) diverges linearly whereas the two Whittaker vec¬ 
tors each vanish linearly, giving a linearly vanishing result in agreement with the behaviour 
of the chiral spinor series (5.36) at s = |. This implies that there is no simplification in 
the Fourier coefficients for s —> 

• In the limit s —>• 2 the Whittaker vector (5.37) remains finite and therefore there is no 
simplification in the Fourier coefficients and the series is of type (2A\)". This is the case 
that is related to adjoint function at s = 1. 

Let us also note that the (2Ai)' type Whittaker vectors always vanish, consistently with the 
analysis of the preceding section. In summary, the chiral spinor series (5.36) always has Fourier 
coefficients attached to the orbits (2Ai)" and A\ (and, of course, the trivial orbit) except for 
s = 1 when the Fourier coefficients reduce to just type A\. 
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5.4 Relation to supersymmetric invariants 

Returning to the V 6 i ? 4 threshold function discussed in this paper and given in (3.1), we conclude 
therefore that the function E|^ooo°] appearing in £(o\) only gets corrections associated to |- 
BPS corrections. Although the adjoint Eisenstein series E^o f o ° J does include constant terms 
inconsistent with perturbative string theory, its wave-front set is of type A 2 , consistent with the 
property that the V 6 i ? 4 threshold function gets corrections associated to -(-BPS instantons. 

The V 4 i ? 4 threshold function is given by a combination of a spinor and vector function 
by [15, 18] 


c( 6) _ 

c (o,i) - 



45 [0 0 0 


(5.40) 


From the analysis of the Whittaker vectors and Fourier modes we now see that the two functions 
have wave-front sets of types ( 2 Ai)' and ( 2 Ai)", respectively, that correspond to the two distinct 
supersymmetric ^-BPS invariants [11]. The general pattern seems to be that ^-BPS corrections 
should be associated with all (special) 2 Ai-type orbits. 
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A Adjoint i?s(8) series as a lattice sum 

In this appendix, we rewrite the adjoint Eisenstein series (2.9) as a sum over a 248-dimensional 
lattice in the adjoint representation. For this we recall that the symmetric tensor product of the 
248 of Eg^g) is 

248 248 = 27 000 0 3875 0 1. (A.l) 

The minimal representation-theoretic constraint to be imposed on a lattice sum is that a charge 
Q E Z 248 C 1R 248 satisfy the constraint that its square contain only the largest representation 
27 000, as is necessary for a lattice sum to be automatically an eigenfunction of the Lapla- 
cian [28]. 

Decomposing the adjoint 248 into 128 © 120 under the Spin( 16) subgroup, one can write 
moduli-dependent es-valued charge VQV -1 (with V the 248-bein) in terms of a Majorana-Weyl 
spinor Xa (with A = 1,..., 128) and the antisymmetric tensor A^ = A^j (with i, j = 1,... , 16) 
of so (16). The invariant mass is then given by XX = XaX a . The representation-theoretic 
constraint that Q <g> Q must have no component in the 1 © 3875 reads 

A ik A jk = ^i(XX) , A ijVX = ipc , Afo-Aju] = ~^{XT ijkl X) . (A.2) 
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By construction the derivative T>a acts on these tensors in the adjoint representation of c 8 ( 8 ), 

i.e. 

V>aX B = ~^AB , DaK] = -TijABX 3 . (A.3) 

Using among other equations that 

4r ijkifi = -4SA^T^-^(Xr^X)T k ^-^(Xr ijM X)-^{XT^ a X)V ijklpqra , (A.4) 
one computes that 

A(AX)" S = V a Va{XX)~ s = 2s(2s - 29)(XX)~ S , (A.5) 

and 


( VT ijkl V)(XX)- s = 2s(2s-5)(XF ljkl X)(XX)- s - 1 , (A.6) 

r kl V(VT ijkl V){XX)- s = 2s(2s-5)(-48(2s-9)A ij + (s-15)T ij f)X(XX)- s - 1 , 

(' Vr i:j T pq V)(VT klpq V)(XX)- s = 8s(2s-5)(96{2s-9){s-7)A ij A kl 

-(s(2s - 41) + 102) (AT*/ kl X) + 2(2s - 21)(s - 22)5 kl {XX)){XX)~ s ~ 1 , 
(VT i[jk r V)(Vr lpq]r V)(XX)- s = -2s(2s - 5)(2s(2s - 29)+48)8 ib {xr klpq] X)(XX )- s ~ 1 . 

We give some indications of how one derives these relations. First, one straightforwardly checks 
that 

Y kl V(VT ljkl V){XX)~ s = s(240o(s)Ajj + 6(s)r^) l(IA)- 8 - 1 , (A.7) 

using representation theory and the 3875 constraint, but computing the explicit coefficients 
a(s ) and b(s) is rather cumbersome. To do so we use the property that the projector to the 
)o i o o o o ° J is defined as 

(: r kl V{VT ijkl V ) + TijV(fA + 336))(XX)~ S = sa(s){ 240A ij - + T ij ^)X{XX)~ s ~ 1 . (A.8) 

This determines b(s). Using moreover the general identity 

T {ij (T™V(VT kl]pq D) + T kl] V(^A + 336)) = -\9(VT ijk {D)(\/A + 180) , (A.9) 

one determines a(s). Once the second identity in (A.6) is known, it is straightforward to compute 
the third one. The last one is more complicated to obtain, but one can nonetheless straightfor¬ 
wardly check that (VT i ] i j k r 'D){VTi pq yV){XX)~ s reduces to the product of a tensor quartic in 
X times (XX)~ S ~ 2 . There is a unique quartic tensor in the )i o o o i o ° j, and a unique quartic 
tensor in the )o o o i o o ° J, and using the 3875 constraint, one obtains that the former vanishes 

(■ XT iyk r X)(XT lpq]r X ) = -8 i{j (XT klpq] X)(XX) , (A. 10) 

such that 

{VT i[jk r V)(VT lpq]r V)(XX)- s = c(s)8 i{j (XT klpq] X){XX)- s - 1 , (A.ll) 


39 


for a coefficient c(s). The latter is fixed using the projector to the |\ o o o 1 o ° J of the fourth order 
differential operator, 


V 


(VT i[jk r V) (: VT lpq]r V ) + <% ( VT klpq] V ) (A + 48) . 


(A.12) 


One finds in this way that the 1/2 BPS equation required for a R 4 type invariant is satisfied 
for s = §, 

(VT ijkl V)(XX)~l = 0, (A.13) 

while the 1/4 BPS equation required for a V 4 i? 4 type invariant is satisfied for s = 

T kl V(VT ijkl V)(XX)-l = -168r^(XX)-f , (A.14) 


and the 1/8 BPS equation required for a V 6 i? 4 type invariant is satisfied for s = -y, 

(VT ibk r V)(VT lpq]r V)(XX)~% = 150 5 i{j (VT klpq] V)(XX)-% . (A.15) 


For any s the function (XX)~ S satisfies that the fourth derivative restricted to the |4 o o o i o ° J 
irreducible representation vanishes. For s = |, the second derivative restricted to the |F o o i o o ° J 
vanishes, for s = | the third derivative restricted to the ^oioooo°J vanishes, and for s = 7 the 
fourth derivative restricted to the [o 2 0 0 0 0 ° j vanishes. 

By analytic continuation, the quartic constraint in the ^1 0 0 0 1 0 ° J is also satisfied by the 
Eisenstein series 


0 l 

[0 0 0 0 0 0 s \ 


1 

2C(2 s) 


E ( X(Q)X(Q))~° , 

Qez 248 


(A.16) 


QxQl38rB=° 


for almost all s. The normalisation here expresses the fact that there is one F/^Z) /-BPS orbit 
on Z 248 for every k G Z. The series only converges for s > 29, but it is defined for complex 
s by Langlands and all the values we consider here are regular. We know that the quadratic 
constraint and the cubic constraint are indeed satisfied by the Eisenstein series at s = | and 
s = |, respectively [17]. However, a more careful analysis exhibits that the Eisenstein series 
does not satisfies the quartic constraint in the |/ 2 0 0 0 0 ° J at s = 7, and the wave-front set does 
not reduce at this specific value. One understands mathematically this property because the 
corresponding nilpotent orbit is not special in this last case. We also note that the adjoint 
series at s = 7 is square integrable and part of the discrete spectrum of the Laplace operator, 
similar to the adjoint SO(5,5) series analysed in detail in section 5.3.2. 


B Poisson summation for the SO( 5, 5) spinor series 

In the string perturbation theory limit, the function E^o 0 0 ° J defined as a sum over pure spinors 
in Z 16 (integers points of an eleven-dimensional variety), decomposes as a sum over chiral and 
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antichiral spinors q and p of Spin(4:, 4) that have a vanishing norm and satisfy (qTp) = 0. This 
gives the sum 


E \ 

0 0 0 


0 

S 



0 ' 
0 


+ 


1 

2 £( 2 s) 


£ £ 

P ez 8 gez 8 

(p,p )=0 (prq)=0 



dt KU\ Z (q +M f +e -<l>\ Z (p)\A 

t l+s e 


(B.l) 


where the lattice sum is divided by 2£(2 s) to get the Langlands normalisation, and = g 
is the six-dimensional effective string coupling constant. One can always find an element of 
Spin(4, 4, Z) to rotate the null spinor p to a preferred basis decomposing as 

8 a = 1 ( “ 2) ® 6 (0) ® 1 (2) , 8 C = 4 (-1) © 4 (1) , 8 = 4 (_1) © 4 (1) , (B.2) 


with respect to SL( 4) C Spin( 4,4), such that 


P = (0,0, rip) , q = (0, < 74 ) , 


(B.3) 


where n p £ Z, and q 4 £ Z 4 . One can therefore use the Poisson summation formula to exchange 
the sum over q with a sum over the dual spinor q = (< 74 , 0 ). Note that the axion decomposes 
as required for the Poisson formula to disentangle the axion dependence as usual, because for 
a = ( 04 , 04 ) one has jp = (0, n p a 4 ) so that q + jp indeed defines a four vector of «SX(4), and 
{p</iq) = gcd(p)(< 74 , 04 ). It remains to compute the determinant for the quadratic term in q. By 
construction, the SL( 4) factor vanishes and one gets the square of the normalised mass square 
of p, ( ^ p ^ ) 2 , such that 
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where we have 


\Z(pTq)\ = \Z(p)\\Z(q)\ . 

At this point it is convenient to use triality, to define the equivalent decomposition 
8 a = 4 (_1) © 4 (1) , 8 C = 4 (-1) © 4 (1) , 8 = 1 (_2) © 6 (0) © 1 (2) , 

such that 


and 


P = ( 0 ,P 4 ) , q = (0, q'4) , (pYq) = (0, 0, (p 4 <?l)) 


|z ® r E g cd(pr«- | '’ w r 2< ‘“ 1) 


|z(p)l‘ 


|^(PP?)| 


(B.4) 

(B.5) 

(B. 6 ) 

(B.7) 

(B. 8 ) 
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where \v(p 4 )\ is the invariant norm associated to the 5L(4) subgroup stabilising (qTp). The sum 
can be replaced by the P^( r L)\SO{ 4,4, Z) coset sum rotating p to the preferred basis (B.2), the 
Pi(Z)\5L(4, Z) coset sum rotating q 4 to a preferred basis and the sum over their respective 
relative integer greatest common divisors. In this basis, (pr<f) = (0, npC/ 4 ), the Pi(Z)\5L(4, Z) 
coset element equivalently determines the direction of ( pTq ) as a vector. It follows by triality 
that one can rewrite this sum as the .Pi(Z)\50(4, 4, Z) coset sum rotating ( pTq ) to the preferred 
basis (B.6) and the Pi(Z)\5L(4, Z) coset sum rotating p 4 to a preferred basis, together with 
the sum over their respective relative integer greatest common divisors, keeping in mind that 
the greatest common divisor of p divides ( pTq ). The Pi(Z)\S'0(4,4, Z) coset sum together with 
the sum over the ( pYq ) greatest common divisor then reduces to the sum over all integral null 
vectors ( pTq ), while 

E 2 E n P 2 KP4)r 2(s-1) = E 2n p {2 ~ s)E [s-io 0 ]{v{ P vq)), (b.9) 

7€Pi(Z)\SX( 4,Z) n p |gcd(prg) n p |gcd(prg) 

where the factor of 2 appears because the sum over n p is then restricted to the positive integers 
dividing ( pTq ). Finally, one obtains after renaming (pTq) as p for simplicity, 


-0 1 — e ^r.nol +e ^ (2s) E \ 00 ° 
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0 0 0 
2 


+ 77TT E (g cd W S_1 E ra2(2 " S) ) 0 0 lK)^-2( 2 7re-^|Z(p)|)e 27ri ^“) (B.10) 

pgZ 8 n|gcd(p) ^ 


(p,p)=0 


where v p parametrizes the SL( 4) C 50(5,5) that stabilises p. 
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